Phase separation in a polarized Fermi gas with spin-orbit coupling 
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We study the phase separation of a spin polarized Fermi gas with spin-orbit coupling near a wide 
Feshbach resonance. As a result of the competition between spin-orbit coupling and population 
imbalance, the phase diagram for a uniform gas develops a rich structure of phase separated states 
involving topologically non-trivial gapless superfluid states. We then demonstrate the phase sepa- 
ration induced by an external trapping potential and discuss the optimal parameter region for the 
experimental observation of the gapless superfluid phases. 
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Spin-orbit coupling (SOC), common in condensed mat- 
ter systems for electrons, has been considered a key in- 
gredient for many interesting phenomena such as topo- 
logical insulators Q], quantum spin Hall effects [2], etc. 
The recent realization of synthetic gauge field and hence 
spin-orbit couplings in ultracold atomic systems opens 
up exciting new routes in the study of these phenomena 
[31 H] , allowing us to take advantage of the features of the 
ultracold atoms, e.g. clean environment and highly con- 
trollable parameters. In particular, with the Feshbach 
resonance technique, the effective interaction strength 
between atoms can be tuned [SJ [3] . This technique has 
been applied to study various interesting topics, e.g. the 
BCS-BEC crossover [7], polarize Fermi gases [8], itinerant 
ferromagnetism [5|, etc. The introduction of spin-orbit 
coupling may shed new light on these strongly correlated 
systems. 

Spin-orbit coupled Fermi gas near a Feshbach reso- 
nance has recently attracted much theoretical attention 
[TUHH)]- The SOC has been shown to enhance pairing 
on the BCS side of the Feshbach resonance [12] , IT4 l [To]. 
Furthermore, for a polarized Fermi gas, the SOC intro- 
duces competition against population imbalance, which 
can lead to topologically non-trivial phases [131 [TH] . Re- 
cently, the phase diagrams for a polarized Fermi gas with 
spin-orbit coupling near a Feshbach resonance have been 
reported for a uniform gas [16] . The phase boundaries 
have been calculated by solving the gap equation and 
the number equations self-consistcntly. However, simi- 
lar to the case of a polarized Fermi gas near Feshbach 
resonance [17] . due to the competition between differ- 
ent phases, the solutions of the gap equation may corre- 
spond to metastable or unstable states. By considering 
the compressibility criterion [16] . the unstable solutions 
are correctly discarded, while the metastable solutions 
may survive, rendering the resulting phase boundaries, 
in particular those representing first order phase transi- 
tions, unreliable. 

In this paper, we examine in detail the zero temper- 
ature phase diagrams for a polarized Fermi gas with 
Rashba spin-orbit coupling near a wide Feshbach reso- 



nance for both the uniform and the trapped cases. To 
avoid getting metastable or unstable solutions, instead 
of solving the gap equation, we minimize the thermody- 
namic potential directly as in Ref. [17] . For the uniform 
gas, we find larger stability regions for the phase sepa- 
rated state at unitarity as compared to the results in Ref. 
|16j . More interestingly, we find that SOC may induce 
more complicated phase separated states involving gap- 
less superfluid phases that are topologically non-trivial, 
in addition to the typical phase separated state composed 
of normal (N) and gapped superfluid (SF) phases. We 
calculate the stability region for the various phase sepa- 
rated states as well as for the gapless superfluid states, SF 
state and normal state. We show that there are two dis- 
tinct gapless phases that differ by the number of crossings 
their excitation spectra have with the zero energy in mo- 
mentum space, consistent with previous results [131 116j . 
These novel gapless phases are stabilized by intermedi- 
ate SOC strengths; whereas for large enough SOC, the 
system always becomes a gapped superfluid of 'rashbons' 
[12] . We show how these phases can be characterized by 
their different excitation spectra and momentum space 
density distributions. We then discuss the phase separa- 
tion in an external trapping potential, where the various 
phases naturally phase separate in real space. By exam- 
ining their respective stability regions, we demonstrate 
the optimal parameter region to observe the gapless su- 
perfluid states in the presence of a trapping potential. 
For all of our calculations in the paper, we adopt the 
BCS-type mean field treatment. Although the mean field 
theory does not give quantitatively accurate results near 
a wide Feshbach resonance, it is a natural first step for 
us to qualitatively estimate what phases may be stable, 
as well as to understand their respective properties. We 
also note that we have neglected the Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) phase in our calculations. This is 
motivated by the fact that the FFLO phase is stable only 
in a narrow parameter region in the absence of SOC due 
to competition against other phases [5j. As SOC intro- 
duces new gapless phases into this competition, we do 
not expect a significant increase in its stability region. 
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We first consider a uniform three dimensional polarized 
Fermi gas with Rashba spin-orbit coupling in the plane 
perpendicular to the quantization axis z. The model 
Hamiltonian takes the form [T31 HH US] 

a k.cr 

k k,k' 

+ ak ±- ( e " llpka l,t a ^-l + h - c ) « (!) 

k 

where £ k = e k — /i, with the kinetic energy e k = 
CT = {til} are the atomic spins; N a denotes the to- 
tal number of particles with spin a; ak j0 -(aj C(T ) annihi- 
lates (creates) a fermion with momentum k and spin tr; 
\i a = /i±/i/2 is the chemical potential of the correspond- 
ing spin species, and V is the quantization volume. The 
Rashba spin-orbit coupling strength a can be tuned via 
parameters of the gauge-field generating lasers [4], while 

fc-L = yj^x + and <^ k = arg (k x + ik y ). In writing 
Hamiltonian ([!]), we assume s-wave contact interaction 
between the two fermion species, with the bare interac- 
tion rate U renormalized following the standard relation 
jj = jj y J^k 2^ H ■ The physical interaction rate is 

given as U p = ^ h m as , where a s is the s-wave scattering 
length between the two fermionic spin species. 

To diagonalize the Hamiltonian, we make the trans- 
formation: a k ,t = ^e i¥,k (a ki+ + a k -), a k 4 = 
^75 («k.+ — flk.-), where a ki ± are the annihilation op- 
erators for the dressed spin states with different helic- 
ities (±) [T2T - HS] . Taking the pairing mean field A = 
V S k ( a -k,iak.t) as m the standard BCS-type theory, 
we may diagonalize the mean field Hamiltonian in the 

basis of the dressed spins: |ak.+ , a^L k + , a k .- : aL k _| • 
The thermodynamic potential is then evaluated from 
n = -ilntr [ e -^ H -^^" N "")], with P = l/k B T. In 
this paper, we will focus on the zero temperature case, 
for which the thermodynamic potential has the form 

n =l E (6-£k,A)-V^, (2) 

k,A=± 



with the quasi-particle excitation spectrum £"k,± = 




Before proceeding, let us examine the quasi-particle 
excitations first and study the conditions for possible 
gapless phases. We see that at the points in the 
momentum space where -Ek,- crosses zero, the quasi- 
particle excitation becomes gapless while the pairing 
gap A remains finite. The SOC, together with the 
population imbalance re-arranges the topology of the 
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FIG. 1. Illustration of typical shapes of the thermodynamic 
potential Q/h as a function of order parameter A/h for var- 
ious phases at unitarity: (a) fi/h = 0.52, ak h /h = 0.1; (b) 
H/h = 0.7, ak h /h = 0.1; (c) fi/h = 0.52, akh/h = 0.3, (d) 
jj,/h = 0.52, akp/h = 0.6. The chemical potential h is taken 
to be the energy unit, while the unit of momentum kh is de- 
ft 2 *; 2 

fined through K = h. 

Fermi surfaces of the spin species [El [16] . The points 
of gapless excitations lie on the k z axis with fcj_ = 0, 
and are symmetric with respect to the k z = plane. 
More specifically, for /i < 0, the excitation spectrum 

has two gapless points ±^pr ( A* + \J *r ~ > 

so long as ^ > yj ' + |A| 2 . For [i > 0, 
the excitation spectrum has four gapless points 

±^(m+ v x-|A| 2 ) 2 ,±^ (/*->/? -|A| 2 ) 2 j 

with |A| < ^ < \J fi 2 + |A| 2 ; two gapless points 

±?f (m + Vt -|A| 2 ) 2 , with f > VM 2 + |A| 2 . We 

identify the superfluid states with two excitation points 
(GP1) and those with four excitation points (GP2) as 
different topological phases [T31 [IB] . 

We illustrate in Fig. [I] typical shapes of the thermo- 
dynamic potential as a function of A with different pa- 
rameters. Notably, due to the competition between dif- 
ferent phases, a double-well structure appears (see Fig. 
[lja-c)). Hence the solutions to the gap equation may 
correspond to the metastable states (local minimum) or 
the unstable states (local maximum) . To make sure that 
the ground state is achieved, we directly minimize the 
thermodynamic potential [17] . 

Another complication comes from the existence of the 
phase separated state, which must be considered explic- 
itly for a uniform gas. As in the case of polarized Fermi 
gases without SOC [8], we introduce the mixing coeffi- 
cient x (0 < x < 1), and the thermodynamic potential 
becomes 

n = xO(Ar) + (1 - x)0(A 2 ), (3) 
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FIG. 2. Zero temperature phase diagram for a uniform Fermi 
gas with population imbalance at (/cfcis) 1 = 0. Within the 
bold phase boundaries are the various phase separated states 
(see text). These phase separated states can be connected 
with the non-phase separated states by first order phase tran- 
sitions (solid bold curve). The thin curves represent various 
second order phase transitions (see text). Here kF = (37r 2 n) 3 , 



E F = 



and n is the total density of the system. 
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FIG. 3. Typical excitation spectrum and momentum space 
density distribution for different phases, (a) Lower branch of 
the excitation spectra for GP1 (solid), GP2 (dashed) and SF 
(dash-dotted) phases; (b-d) Density distribution in momen- 
tum space for spin-up (solid) and spin-down (dashed) species 
along k± — and k z = (inset), for (b) akh/h = 0.35, 
fi/h = 0.5 (SF); (c) ak h /h = 0.35, fi/h = 0.45 (GP1); (d) 
ak h /h = 0.45, fi/h = 0.43 (GP2), respectively. 



where (i = 1,2) is the pairing gap for the ith com- 
ponent state. Note that due to SOC, we now have the 
possibility of a phase separated state of two distinct su- 
perfluid states (see Fig. [ljc)). The number equations of 
the phase separated state become 



(1 



A=Ai 



(4) 



A=A 2 



Minimizing the thermodynamic potential Eq. (|3| with 
respect to A, and x while implementing the number con- 
straints Eq. Q , we map out the phase diagram for a uni- 
form polarized Fermi gas with SOC at (kpas)^ 1 = 0. Fig. 
[2]illustrates the resulting phase boundaries in the plane of 
(P, akp/Ep), where the polarization P = jf-^f ■ When 
the SOC is off (a = 0), the system remains in a phase 
separated state of normal and gapped superfluid (PS1) 
up to P ~ 0.93 before it becomes a normal state via a 
first order phase transition. This is consistent with pre- 
vious mean field calculations for a polarized Fermi gas 
[51 HH] , while different from the result in Ref . [T5] . As 
the SOC strength a increases, a rich structure of dif- 
ferent phases shows up, e.g. gapped superfluid phase 
(SF), gapless superfluid phases with different Fermi sur- 
face topology (GP1 and GP2), and notably, various phase 
separated states. These phase separated states are con- 
fined by a phase boundary of first order phase transition 
(bold curve in Fig. [2J). In addition to the typical PS1 
phase, we now have a phase separated state with GP2 
and SF phases (PS2), and a phase separated state of two 
distinct GP2 phases (PS3). As a increases, the system 
can undergo second order phase transitions from PS1 to 
PS2 and then to PS3 for intermediate P and a. As- 
suming | Ai| < IA2I, the phase boundaries between them 



can be determined by imposing Ai = (PS1 and PS2) 
and I = |A 2 | (PS2 and PS3), respectively. These phase 
separated states finally become unstable and give way to 
single component superfluid phases as a becomes large. 
The phase boundaries between these single component 
states are determined by setting | = |A| (SF and GP2), 
I = V-" 2 + |A| 2 (GP2 and GP1), and A = (N and 
GP1), respectively. When a is large enough, the stabil- 
ity region of the GP2 phase decreases and finally vanishes 
at a tri-critical point (/i = 0), beyond which only GP1, 
SF and normal phase may exist. Note that beyond the 
tri-critical point, the chemical potential fj, becomes neg- 
ative, and the phase boundary between GP1 and SF will 
bend upwards so that in the large a limit the SF phase 
becomes dominant in the phase diagram. 

To characterize the properties of the different 
phases, we calculate the excitation spectrum and 
number distribution in momentum space for SF, 
GP1 and GP2 states (see Fig. |3|. Several interest- 
ing observations are in order. Firstly, the gapless 
phases leave their signatures in the momentum 
space density distribution. For fcj^ = and \k z \ G 




the minority spin component vanishes, and pairing does 
not occur in this region. This is reminiscent of the 
momentum space phase separation of a breached pairing 
phase in the polarized Fermi gas [TS], though now the 
unpaired region lies only on the k z axis. Away from 
k z axis, the occupation of the minority spin recovers 
from zero gradually, leaving a signature which may be 
detected in the time of flight imaging experiment [16] . 
Secondly, for finite a, both the gapless and the gapped 
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FIG. 4. Phase diagram in the (p/h, aku/h) plane at 
(fcfta s ) _1 = 0. While the second order phase transitions are 
in dashed thin curves, the first order phase transitions are 
shown in solid bold curves, which end at the point where the 
double-well structure in the thermodynamic potential disap- 
pears (inset). The boundary for vacuum (V) is determined 
by setting the chemical potential of the majority spin species 
to vanish in the normal phase. 



superfluid phases can support population imbalance, 
which can be seen from the density distribution along 
k± (see Fig. [3] insets) . Indeed as we will see later, for 
large enough a, we may expect no phase separation even 
in the presence of a harmonic trapping potential. The 
atoms in the trap will all be in the superfluid phase 
induced by SOC. 

To understand the spatial distribution of the various 
phases in a trapping potential, we calculate the phase 
diagram as a function of (p/h, ak^/h) at unitarity (Fig. 
[I]), where kh is defined in the caption of Fig. [I] Under 
the Local Density Approximation (LDA) while assuming 
both spin species experience the same harmonic poten- 
tial, the local chemical potential jit(r) can be related to 
that at the center of the trap [i as /x(r) = fi — V(r), 
where V(r) gives the trapping potential. Thus a down- 
ward vertical line in Fig. [4] represents a trajectory from 
a trap center to its edge, with the chemical potential at 
the trap center fixed by that at the starting point of the 
line. In Fig. |4j consistent with Fig. [2j the GP2 phase 
only exists in a small parameter region in the trap, while 
there appears to be considerable stability regions for the 
GP1 phase. When a is small, the Fermi gas in the trap 
will phase separate into two regions, SF at the core, nor- 
mal phase (N) towards the edge. At intermediate a, the 
gapless phases GP2 and GP1 may appear either near the 
center of the trap or as a ring between the SF core and 
the normal edge, depending on the chemical potentials. 
Note that the boundary of the first order phase transi- 
tion between PS3 and GP2 (dotted thin curve in Fig. [2]) 
corresponds to a small scale structure here (Fig. |4]insct), 
where a first order phase transition (bold black curve) ex- 
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FIG. 5. Phase diagram in the (^j,/h,akh/h) plane at (a) 
(khCis)^ 1 = —1 and (b) (fe^eis) -1 = 0.5. First order phase 
transitions are shown in solid bold curves, while second order 
phase transitions are in dashed thin curves. 



ists between two distinct gapless superfluids, both in the 
GP2 phase. However, this region is found to be small at 
unitarity and only increases slightly towards the BCS side 
of the resonance. It is therefore difficult to observe this 
phase transition in the trapped Fermi gas in the parame- 
ter region that we considered. For large SOC beyond the 
tip of the GP1-SF phase boundary, there is only the SF 
phase in the phase diagram, and hence we will have only 
the SOC induced SF phase in the trap for large enough 
SOC. 

We have also calculated the phase diagram in 
{l-i/h, akh/h) plane away from the resonance point. On 
the BCS side (Fig. ga)), the stability region for the GP2 
phase increases considerably. It is therefore desirable to 
prepare the system on the BCS side of the resonance for 
the observation of GP2 phases. On the BEC side (Fig. 
[5^b)), the GP2 phase vanishes from the phase diagram 
altogether for \i < 0, consistent with our previous discus- 
sion. 

In summary, we have calculated in detail the phase di- 
agrams near a wide Feshbach resonance for a polarized 
Fermi gas with Rashba spin-orbit coupling. We find that 
the competition among pairing, polarization and SOC 
gives rise to a rich structure of phases and phase separa- 
tions involving topologically non-trivial phases. From the 
phase diagrams for both uniform and trapped systems, 
we find that the interesting gapless superfluid phases are 
most likely to be observed in an experiment with moder- 
ate polarization and SOC strength. 
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